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ABSTRACT 

We have used the WFPC2 camera on the Hubble Space Telescope to obtain 
photometry of the outer-halo globular clusters Palomar 3, Palo mar 4, and 
Eridanus. These three are classic examples of the "second parameter" anomaly 
because of their red horizontal-branch morphologies in combination with their 
low-to-intermediate metallicities. Our color-magnitude diagrams in (V, V-I) 
reach Vu m ~ 27.0, clearly delineating the subgiant and turnoff regions and about 
three magnitudes of the unevolved main sequences. The slopes and dereddened 
colors of the giant branches are consistent with published [Fe/H] estimates that 
rank the clusters (Pal 3, Eridanus, Pal 4) in order of increasing metallicity, with 
all three falling near or between the abundance values of the classic nearby halo 
clusters M3 and M5. Differential fits of their color- magnitude diagrams are 
made to each other and to M 3 and M 5 for relative age determinations. We find 
that the three outer-halo cluster CMDs differ from the nearby clusters in a way 
that is consistent with their being younger by ~ 1.5 — 2 Gyr, if we have correctly 
estimated the clusters' chemical-abundance ratios. Conversely, the inferred age 
difference could be smaller ( < 1 Gyr) if either [Fe/H] or [a/Fe] for the outer-halo 
clusters is significantly lower than we have assumed. Possible age spreads of 
order 1 Gyr among both the nearby and outer-halo clusters may also be present. 

Subject headings: Globular clusters: individual (Palomar 3, Palomar 4, 
Eridanus); Galaxy: formation 
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1. Introduction 

This is the second paper in which we report results of age measurements for the 
globular clusters at the outermost reaches of the Milky Way halo. In the first paper of this 
series (Harris et al. 1997 [Paper I]), we showed that the massive globular cluster NGC2419, 
at a Galactocentric distance of ~ 90 kpc, has a color-magnitude diagram ("CMD") that is 
indistinguishable from that of M 92 and other nearby halo clusters with the lowest known 
metallicities, [Fe/H] ~ —2.2. These clusters have the same age to within the measurement 
precision of ~ 1 Gyr, suggesting that globular cluster formation may have started at 
virtually the same time everywhere across the entire expanse of the protoGalaxy. 

However, NGC2419 represents only the upper envelope of the chronology of the 
outermost-halo clusters. The other remote clusters (Palomar 3, 4, 14, Eridanus, and AM-1) 
are intriguing in quite a different way as strong cases of the classical "second parameter" 
problem: their horizontal-branch ("HB") morphology is redder than the mean for nearby 
halo clusters of similar (low to intermediate) [Fe/H]. The tendency for the second-parameter 
effect to become increasingly evident in the outer halo is one of the key motivations behind 
the Searle & Zinn (1978 [SZ]) scenario in which the Galaxy's spherical halo is envisaged 
as having accumulated over time in a series of mergers of sub-galactic fragments. Making 
the plausible assumption that age is the most important second parameter responsible 
for intrinsic scatter in the relation between [Fe/H] and HB morphology, SZ concluded 
that the mean age of the clusters decreases, and the age spread increases with increasing 
Galactocentric distance, -Rgc- By contrast, if all parts of the spherical halo formed into 
stars during a rapid free-fall collapse as was argued by Eggen, Lynden-Bell & Sandage (1962 
[ELS]), it would be much more likely that the halo clusters would have rather similar ages 
and would (for example) follow a simple relation between CMD morphology and metallicity. 

Arguing from a comparison of quantitative simulations of HB morphology as a function 
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of metal abundance and age to CMDs of real clusters Lee, Demarque & Zinn (1994, 
hereinafter LDZ; see especially their Fig. 7) concluded that the globular clusters within 
~ 8kpc of the Galactic center — spanning the full metal-abundance range from less 
than one one-hundredth solar to nearly the solar value — are essentially coeval. Those at 
Galactocentric distances between 8 and 40 kpc formed an average of some 2 Gyr later, but 
with a perceptible dispersion in their formation epochs. Finally, the five Galactic globulars 
with red horizontal branches that lie outside 40 kpc formed more than 4 Gyr later than their 
counterparts in the Galactic center (these statements are based upon our interpretation of 
LDZ's Fig. 7). 

Lee (1992) has added another wrinkle to the earliest formation history of the Galaxy, 
using another argument based on HB morphology. He notes that the RR Lyrae variables 
in Baade's window — a region where we are allowed a more or less clear glimpse into the 
Galactic bulge — possess metal abundances closely bunched about a value of [Fe/H] = -1.0. 
He then argued that for stars of such high metal abundance to populate the HB instability 
strip, they must be ~ 1 Gyr older than even the mean age of the globular clusters with 
i?cc < 8 kpc. 

However, the foregoing arguments rely on the assumptions that age is the dominant 
second parameter after metal abundance, and that HB morphology can be used in a simple 
and direct way as a measure of cluster age. These assumptions are straightforward and 
widely accepted, but they are perhaps not yet ironclad, and the debate continues (e.g., 
Stetson, VandenBerg & Bolte 1996; Sarajedini, Chaboyer & Demarque 1997). There are 
two additional steps required to test and calibrate this chronology directly: first, we must 
measure in considerable detail the chemical compositions of stars in these clusters (since 
composition differences can mimic age differences in a variety of ways); and second, we 
must make much more direct estimates of age based on photometry of the main sequence 
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and subgiant branch in clusters spanning the range of Galactic environments. Stars in these 
stages of evolution are driven by somewhat simpler internal physics, and suffer from fewer 
potential complications due to uncertain mid-life processes such as mixing and mass-loss, 
than the HB stars (see, e.g., Sweigart 1997). 

Among the six Galactic globular clusters at Galactocentric distances greater than (to 
choose an arbitrary boundary) that of the Large Magellanic Cloud, only one — Pal 14 - 
has had its main-sequence turnoff detected from the ground (Sarajedini 1997), although 
even in this best case the data penetrate less than a magnitude past the turnoff. Sarajedini 
argues that Pal 14, with [Fe/H] -1.6, is 2 Gyr younger than the nearby globular M5 
and 4 Gyr younger than NGC 6752. Additional caveats on these results are that M5 may 
have a significantly higher [Fe/H] than Pal 14, and NGC 6752 has an extremely blue HB 
morphology that makes it very difficult to register it and Pal 14 to a common distance 
modulus. 

It is possible that the outermost-halo objects may actually bear little direct relevance 
to the formation history of the main part of the Galaxy. Many or all of the red-HB clusters 
in particular (Pal 3, 4, 14, Eridanus, AM-1) may belong to a stream of remote satellites 
(the Fornax-Leo-Sculptor stream; see Majewski 1994) in an orbital plane with an extremely 
large semimajor axis, R ~ 150 kpc or more. No matter whether the Galaxy itself formed in 
a monolithic or a hierarchical fashion, it is possible that this handful of objects was born 
well apart from the original protoGalaxy, and are now left orbiting a large stellar system of 
which they were never really an integral part. If this is true, these clusters may merely be 
a red herring in the search for clues to the formation of the Milky Way. But at least these 
same clusters should still provide us with a test of the proposed correlation between HB 
morphology and age, which is currently regarded as the essence of the second-parameter 
problem. In this paper, we use data from the Hubble Space Telescope to address the 
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questions of whether, and to what extent, the intermediate-metallicity clusters of the outer 
halo and the nearby halo really do have different ages. (Following LDZ, we separate the 
"inner" halo from the "intermediate" halo at a Galactocentric distance of about R & , and 
the "intermediate" from the "outer" halo at about 50 kpc. In this paper, we use the word 
"nearby" in its common English sense, to describe those clusters that are close to the Sun 
and, hence, usually quite well studied. It encompasses some of the outermost members of 
the inner halo as well as the closer members of the intermediate halo; "nearby clusters" 
without the additional modifier "halo" can also include a few of the nearer disk and/or 
bulge clusters.) 



2. Observations and Data Reduction 

The clusters we report on here are Palomar 3, Palomar 4,| and Eridanus. These are 
three of the six known globulars beyond a Galactocentric distance Rqc ~ 50 kpc, and 
(typically for such objects) have low luminosities, diffuse structures, and low to intermediate 
[Fe/H] values. Pal 3 and Pal 4 were imaged with the WFPC2 camera during HST Cycle 
4, with Pal 4 being visited twice. Eridanus was observed during Cycle 5. As was the case 
for NGC2419, the first visit to Pal 4 resulted in eight long exposures in F555W (each 
of 1400 s duration) and F814W (of 1300 s), after which the spacecraft was shifted by 20" 
toward the center of the cluster, and six short exposures were taken in each filter (60 s 
in F555W and 30 s in F814W). For the second visit to Pal 4 and the visit to Pal 3, a 
different scheme was employed. For these observations we specified that the cluster be 

2 Results from reliminary analyses of the present data for Pal 3 and Pal 4 were included in 
the study by Richer et al. (1996). The current analysis supersedes the results cited in that 
paper. 
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centered on the PCI chip and requested cycles of exposures, each cycle consisting of an 
(approximately) 1800 s, a 60 s, and a 30 s exposure; we obtained eight such cycles per filter 
for each cluster. Because Eridanus is closer than the other two clusters, for this system 
we obtained seven exposures of 1100 or 1200 s duration in each of F555W and F814W, 
plus 3 x 180 s in F555W and 3 x 230 s in F814W. Extensive additional details about these 
observations may be obtained by electronic query to either the Space Telescope Science 
Institute in Baltimore ( |http: / / archive. stsci.edu/HDA /| ), or the Canadian Astronomy Data 



Centre ( http: / / cadcwww.hia.nrc.ca/hst / scicncc.html ; specify GO programs 5481, 5672, and 



6106: J. E. Hesser, PI), or by contacting the first author of this paper. 

The process of data reduction was the same as that for NGC 2419 in Paper I, and will be 
described more completely in Stetson (1998, in preparation). The suite of codes collectively 
referred to as the third-generation DAOPHOT (including ALLFRAME; see Stetson 1994) 
was employed to obtain instrumental magnitudes and colors for all measurable stars in 
the WFPC2 fields and to merge these onto a common photometric system. Raw observed 
instrumental magnitudes were corrected for charge loss due to charge-transfer inefficiency 
by an amount equal to -0.04 mag per 800 pixels times the ^-coordinate of the star on the 
chip (Holtzman et al. 1995; Whitmore & Heyer 1997). Color terms for the transformation 
from the instrumental (F555W, F814W) magnitude system to the standard Johnson V, 
Kron-Cousins / systems were taken directly from Holtzman et al. (their Table 7). The 
zero-point calibration relies ultimately on ground-based images of some of these same fields 
obtained by us. We found that the zero points inferred for our short exposures agree with 
those from Holtzman et al. to within ±0.02 mag in both V and /. However, the zero points 
appropriate to our long exposures appear to differ from those of Holtzman et al. by some 
0.05 mag in each filter, in the sense that the effective quantum efficiency of the CCDs in 
WFPC2 appears to be higher in long exposures than in short ones. We believe that this 
behavior is consistent with the notion that crystal-lattice imperfections in the CCDs act 
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as charge traps. In long exposures, where the diffuse sky foreground is appreciable, these 
traps consume electrons and lower the overall diffuse flux in the image without affecting 
stellar profiles or brightnesses. However, in short exposures the diffuse sky is not sufficient 
to fill all the traps, and electrons are consumed from the charge packets representing 
stellar images either during the exposure itself or during readout of the chip. This effect 
appears to operate in addition to the charge-transfer inefficiency ramp of -0.04 mag per 
800 pixels mentioned above. Fortunately, for the differential comparisons to be performed 
below accurate knowledge of the zero points is not essential: incorrect values will slightly 
falsify the distances and/or reddening values inferred for our clusters, but the relative age 
indicators will be unaffected. We hope to be able to improve the fundamental accuracy of 
our zero points once the charge-transfer inefficiency of the WFPC2 CCDs becomes better 
characterized. For our present purposes, separate photometric zero-points for the long 
and short exposure times have been determined by direct comparison to our ground-based 
observations of Pal 4 and NGC2419. 

Later, we will compare the results for these three clusters to our ground-based 
photometry for the nearby halo clusters M 3 and M 5. Data for M 3 were obtained during 
seven observing runs on four different telescopes between 1983 and 1994 (although the 
/-band data came from only two of those runs), while M5 was observed during twelve runs 
on six different telescopes during the same time span (again the /-band data were from two 
of those runs). These data include the observing run analysed by Johnson & Bolte (1998), 
but the images have been independently reanalysed and recalibrated for the present study. 
All the data from the various runs were homogenized to a common photometric system - 
that of Landolt (1992) — as originally outlined some years ago by Stetson (1990, 1993). 
The present ground-based photometric system is anchored to a total of some 44 observing 
runs in V and 17 runs in /. 
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3. The Color-Magnitude Diagrams 

Figs. 1-3 Figures 1 through 3 show the CMDs we have derived from the HST imagery for 
the remote clusters, Pal 3, Pal 4, and Eridanus. For each object, the data from the four 
WFPC2 chips have been combined. The narrowness of the red giant branches ("RGBs") 
in all three CMDs is consistent with our expectation that the four WFPC2 CCDs and 
the long and short exposures have all been normalized to the same photometric system to 
within ~ ±0.01 mag. We note that all three of the outer-halo clusters contain significant 
numbers of blue straggler stars (about a dozen such objects appear in each of Figs. 1-3). 
Figs. 4 & 5 Figures 4 and 5 illustrate the ground-based color- magnitude diagrams for M 3 and M 5 that 
we will later compare to those of the outer-halo clusters. The present data for these clusters 
give results that are very similar to those of Johnson & Bolte (1998), with about twice the 
sample size. 

Palomar 3 is the only outer-halo cluster in our study which appears to contain RR 
Lyrae variables. Three were found by Gratton & Ortolani (1984), who also identified one 
candidate Population II Cepheid; a fifth variable candidate, which had previously been 
noted by Burbidge & Sandage (1958), was considered to be non- variable by Gratton and 
Ortolani. One of Gratton & Ortolani's RR Lyrae stars, their number 34, lies outside our 
field coverage, as does the candidate Cepheid (their number 102). We have recovered the 
other two candidate RR Lyraes from Gratton and Ortolani; their number 155 = our number 
1-233, and their number 283 = our number 2-614. In addition, in contrast to the conclusion 
of Gratton and Ortolani, we find that the Burbidge-Sandage candidate = Gratton- Ortolani 
number 188 = our number 2-198, is indeed a real variable. Finally, we identify four 
Fig. 6 additional RR Lyrae candidates not previously noted. In Figure 6 we present light-curve 
fragments from our data for these seven stars. Within the limitations of the available data, 
all seven of the variables appear to have the asymmetric sawtooth lightcurve characteristic 
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Table i of RR a 5-type variables. Table 1 represents notional period estimates for the stars from 

these data (for our present crude purposes, we neglect the slight difference in phase between 
extrema as observed in the two photometric bandpasses). We think it most likely that for 
five of the seven candidate variables, three full cycles occurred between the two dates of 
observation, while for the other two, candidates 1-299 and 2-614, barely more than two 
cycles elapsed between the end of the first observing sequence and the start of the second. 
These assumptions yield virtually identical periods of 0.6 d for all seven of the RR Lyrae 
candidates. The possibility of four cycles having occurred for any of the variables, implying 
a period of order 0.44 d, can be ruled out by the absence of significant lightcurve overlap on 
either occasion, since on each date we obtained continuous data strings more than 0.48 d 
in length. Conversely, the possibility that some of the variable candidates underwent only 
two full cycles or less during the course of observations would imply periods of order 0.9 d, 
which are extremely rare among classical RR Lyraes. The mean magnitude and color of 
all seven candidates is indicated in Fig. 1 as (V) = 20.49 ± 0.027, ( V-I) = 0.66 ± 0.061. 
Among the non-variable stars in our Pal 3 sample, we count six in the neighborhood of the 
red horizontal branch, implying a horizontal-branch morphology index for this cluster of 
(B - R)/{B + V + R) = -6/13 = -0.5 ± 0.2. 

We have determined fiducial points for the principal CMD sequences in two ways. 
First, we simply calculated robust mean magnitudes and colors (after three iterations of 
sigma-clipping) in bin steps of Ay ~ 0.2 mag for the three outer-halo clusters, and in steps 
of 0.1 mag for M3 and M5. These mean points are shown as squares in Figures 7 through 
11 and are listed in Table 2: columns 1 and 2 give the mean (V, V — I) values for the stars 
found in each bin, while the last three columns give the number of stars per bin retained 
in each of the three iterations. Notice that the bin steps are not strictly uniform in AV, 
particularly along the upper giant branch, where the bin intervals and sizes are determined 
more strongly by the (patchy) distribution of RGB stars. 
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Second, we took D. A. VandenBerg's latest set of theoretical isochrones and, regarding 
them as a set of numerical drafting splines likely to resemble an actual cluster locus more 
closely than, say, a low-order polynomial, fitted them directly to the original data. These 
isochrones are essentially the ones described in Paper I except that they now incorporate the 
latest neutrino cooling rates (Itoh et al. 1996), which result in a small increase in the derived 
core mass at the giant-branch tip (0.002-0.004 M ) and an increase of ~ 0.025 mag in the 
corresponding HB luminosity. These changes are irrelevant for our immediate purposes. 
Some complications, such as the uncertain effects of convection and mixing on the derived 
radii for evolved stars, reddening, and the practical difficulty of accurately characterizing 
the relationship between colors on the Johnson/Kron-Cousins standard system and those 
obtained with a particular filter /detector combination, continue to constitute problems 
for the reliable absolute interpretation of CMDs (e.g., VandenBerg 1998, in preparation). 
The magnitude and color zero-points of the isochrones employed here are based on the use 
of theoretical zero-age HB models as absolute luminosity fiducial points, but it turns out 
that our models are numerically consistent with a compromise between cluster distance 
estimates from the latest Hipparcos parallaxes for subdwarfs on the one hand (e.g., Pont et 
al. 1998), and for RR Lyraes on the other (Fernley et al. 1998). Clearly, absolute cluster 
ages will be quite a strong function of one's assumptions concerning the distance scale (see, 
e.g., Chaboyer et al. 1998 and Gratton et al. 1997; also, Hendry 1997), but unless we are 
very unlucky age differences from direct cluster-to-cluster comparisons will be minimally 
affected. 

We used theoretical isochrones for [Fe/H] = -1.14, -1.31, -1.41, -1.54, and -1.61; 
[a/Fe] = 0.0, +0.3, and +0.6; and age = 8, 10, 12, 14, 16, and 18 Gyr. We numerically 
fitted all of these isochrones to the data for each of the clusters, allowing arbitrary vertical 
and horizontal shifts in each case. Isochrones were fitted to the actual cluster data (i.e., not 
to the normal points), and the fits were carried out in the (I, V-I) plane, because in this 
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representation the subgiant branch is more strongly sloped than in the (V, V-I) diagram. 
Looked at in another way, if apparent magnitude is taken as the independent variable and 
color as the dependent one, then the long, flat subgiant branch characteristic of young ages 
becomes vertical in the function V-I = Fn(V); indeed, for very young ages the relationship 
can be multi-valued so the mapping between color and magnitude becomes equivocal. The 
relation V-I = Fn(7) is better behaved in this regard. 

The fits were carried out using a robust maximum-likelihood scheme that considered 
the photometric uncertainty of each stellar measurement and reduced the relative weights 
of obvious outliers. This method (essentially the same as the one illustrated in Fig. 6 
of Stetson 1987) is quite insensitive to the actual numbers and distribution of outliers, 
and is also robust against variations in the details of the weighting scheme. After all the 
isochrone fits had been done, the one isochrone which produced the highest maximum of 
this likelihood function was identified and adopted as representing the position and shape 
of the turnoff and subgiant region of the CMD. Note that in this procedure we imposed 
no a priori constraints on cluster distance, reddening, chemical abundance, or age within 
the parameter space spanned by these particular isochrones, nor do we intend to draw 
any immediate a posteriori conclusions concerning those parameters. At this point, all 
we want is to utilize the theoretical isochrones as a set of numerical curves describing the 
shape and position of the turnoff and subgiant branch in the CMD. That said, we note 
that the specific isochrones selected by the statistical procedure were, in fact, reasonable: 
the best matches were achieved with model metal abundances in the range -1.41 < [Fe/H] 
< -1.61 and model ages ~ 14-16 Gyr. We stress again that no great weight is to be placed 
on these numerical quantitities, apart from a general sense of satisfaction that they are not 
unreasonable. 

Having determined which isochrone best matches the appearance of the cluster 
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sequence within two magnitudes of the turnoff, Ito — 2 < / < I TO + 2 mag (solid curves in 
Figs. 7-11), we then read out from the isochrones the magnitudes and colors corresponding 
to five points: (1) the turnoff, (2) the point on the subgiant branch that is precisely 
0.10 mag redder than the turnoff, (3) the point on the subgiant branch that is 0.05 mag 
redder than the turnoff, and (4) and (5) the points on the upper main sequence that 
are precisely 0.05 and 0.10 mag redder than the turnoff. After the (/, V-I) locations of 
these fiducial points had been determined, the magnitudes were converted to their V-band 
equivalents. We present the apparent visual magnitudes of isochrone points (2) through 
Table 3 (5), denoted by V^7 010 , V+o.o5> V+omi an d V+o.ioi respectively, in Table 3. For purposes of 
the present discussion, we took the turnoff of an isochrone to be given by the magnitude 
Mi jT o such that the (V-I) color of the isochrone at a point 0.2 mag more luminous than 
M I: to is exactly equal to the color at a point 0.2 mag less luminous than M I:T o- The 
tabulated color on the fitted isochrone at Mj^o was then taken as the turnoff color of the 
cluster. This definition was adopted to avoid uncertainty due to minor numerical glitches 
either in the theoretical evolution calculations or in the interpolation from evolutionary 
tracks to isochrones. However, we subsequently found that the turnoff color defined in this 
way agreed to within 10~ 4 mag of the bluest color actually achieved on the isochrone, so it 
seems our concerns were exaggerated. 

VandenBerg, Bolte & Stetson (1990, "VBS") suggested using the main-sequence point 
"V+o.05 to effect the vertical registration of cluster sequences for estimating relative ages, 
while Chaboyer et al. (1996a) advocated the use of the subgiant-branch point V+ 0m05 . 
The advantages of estimating a fiducial point on the subgiant branch are, first, that 
the sequence is more nearly horizontal than the main sequence, and second, that the 
photometry for individual stars is likely to be more precise. A fiducial point on the main 
sequence has the advantage that it can be estimated from a much larger number of stars. 
Furthermore, the shape of the upper main sequence is "simpler," both in the sense that 
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the high-order derivatives are smaller and in the sense that the overall shape of the upper 
main sequence is less sensitive to age, abundance, and uncertain details of stellar interiors 
than the morphology of the subgiant branch (e.g., VBS Figs. 2 and 3). However valid these 
distinctions may be in general, with the present data sets both the subgiant-branch and 
main-sequence fiducial points lead to much the same relative shifts, as Table 3 shows, and 
hence ultimately to the same astrophysical conclusions. In fact, since the whole shape of 
the turnoff region has been fitted as a fixed unit, when the various fiducial points have been 
determined in this way all are determined to essentially the same level of precision, and it 
makes little difference whether subsequent analysis depends upon the turnoff magnitude or 
upon either of the subiant-branch or main-sequence fiducial points. 

We note that the isochrone fits imply turnoff colors approximately 0.007 mag bluer than 
the bluest normal points given in Table 2; possibly this is because the nose of the turnoff is 
more pointed in the models than in reality, but it is also possible that it is a consequence of 
a different level of effectiveness in rejecting binaries and uncertain data in the two schemes. 
In either case, since it appears to be a constant effect (at a level of ±0.003 mag or so), it 
is unimportant for differential comparisons among clusters provided care is taken to base 
those comparisons on one type of analysis at a time. Apart from this minor difference, 
we see generally excellent agreement between the traditionally defined normal points and 
the somewhat more novel isochrone fits, indicating that the latter method does provide a 
robust, impersonal, and repeatable method for establishing the position and morphology of 
a cluster locus from unbinned data. 

3.1. Intercomparison of the Outer Halo Clusters 

The first-order variations of our present isochrones with changes of age, overall metal 
Figs. 12-14 abundance, and a-element enhancement are illustrated in Figures 12-14, where we have 
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normalized a selection of isochrones according to the precepts of VBS: the isochrones are 
translated horizontally to match their turnoff colors (vertical short-dashed line), and then 
they are shifted vertically to register the point on the upper main sequence that lies exactly 
0.05 mag redder than the turnoff (cross). The variation of isochrone morphology with a 
change in one of the assumed input parameters is then most easily perceived as a change in 
the position of the giant branch. For subsequent comparison with our cluster data, we have 
parameterized these morphology variations by the color difference between the turnoff and 
the point on the giant branch (horizontal long-dashed line) that is precisely 3.0 mag brighter 
in V than the upper-main-sequence fiducial registration point (cross); this corresponds to 
a magnitude level 2.19 or 2.20 mag brighter than the turnoff, as we have defined it above. 
We will represent this color difference by (V-I)gb — (V-I)to, or A(V-I). Clearly, when 
isochrones are compared to cluster data or cluster sequences are compared to one another 
after registration in this fashion, the differential interpretations that one may then make 
are independent of cluster distance and reddening and of zero-point errors in both the 
photometry and the theoretical isochrones. We find that the change in giant-branch color 
associated with a change in presumed age r from 8 to 16 Gyr is 5[A(V-I)] = —0.1 15 mag, 
or 5[A(V-I)]/5\ogT = —0.382 mag/dex, which — in first-order expansion about an age of 
12 Gyr — corresponds to about -0.014 mag/Gyr. Similarly, to first order A(V-I) decreases 
by about 0.04 mag for a 1 dex increase in either [Fe/H] or [ce/Fe]. 

The true cluster abundances [Fe/H] and [a/Fe] are therefore significant for interpreting 
these data, and we discuss current estimates in detail in the Appendix. The literature 
values are not definitive, but they suggest that the three clusters are quite similar in 
[Fe/H], but with Pal 3 likely to be the most metal-poor, Pal 4 the most metal-rich, and 
Eridanus nearly as metal-rich as (and perhaps indistinguishable from) Pal 4. In Figure 15 
Fig. 15 we compare the normal points of the three clusters after dereddening and registration 
of their HB levels (details to be discussed below). The reddening values E(B-V) are 
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taken from Harris (1996), "Catalogue of Milky Way Globular Cluster Parameters" (see 
[tittp: / / www.physics.mcmaster.ca/Globular.html| , revision of 1997 May 15), and in our 
present analysis we have taken E(V-I) = 1.3E(B-V). (This value comes from interpolation 
within Table 3 of Cardelli, Clayton & Mathis 1989, taking V ~ 555 nm and I ~ 814 nm; 
recall that the reddening actually occurs in the instrumental magnitudes and not in the 
standard ones. The reddening values estimated for these five clusters are so small and 
so similar that our results are quite insensitive to the value assumed for this ratio.) For 
purposes of this illustration only, we have assumed that all three clusters have an HB 
luminosity My = +0.70 at (V-I) = 0.60 (a color close to that of the turnoffs). This 
photometric comparison supports the conclusion that the clusters are indeed similar in 
[Fe/H] and that the relative metallicity ranking from the literature is probably correct: 
Pal 3 has the steepest and bluest RGB, while Eridanus and Pal 4 have RGBs with slightly 
redder colors and flatter slopes that are virtually identical to one another. In addition, the 
HB of Pal 3 extends blueward into the RR Lyrae instability strip, whereas the HBs for both 
Eridanus and Pal 4 are simply red stubs with extremely small ranges in color, which — all 
other things being equal — also tends to corroborate the inferred metallicity ranking. 

The near-coincidence of the main-sequence turnoff ("MSTO") and subgiant regions 
Fig. 16 among all three clusters immediately suggests that they have similar ages. Figure 16 shows 
the three cluster fiducials after registration via the VBS method. After such registration, 
any age differences should be manifest as offsets in the cluster RGBs (again, all other things 
- such as composition — being equal). None is obvious in this diagram. However, we can 
Fig. 17 perform a more quantitative test: Figure 17 shows the color-magnitude diagrams of Pal 3, 
Pal 4, and Eridanus with the three clusters registered to match their apparent turnoff colors 
and fiducial magnitudes derived from the isochrone fits (use of the fitted turnoff 

magnitudes or the V+ 05 fiducial points in effecting the vertical registration would have 
made no material difference to the comparison). To eliminate some poorer measurements, 
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here we consider only those stars with photometric uncertainties a(V-I) < 0.05 mag, 
where the standard error assigned to a given star represents a compromise between that 
derived from the expected Poisson and readout noise and that indicated by the actual 
frame-to-frame repeatability for that particular star; 13% of the measured stars were 
rejected by this criterion. After registration of the data, we fit a simple parabola to the 
color of the giant-branch stars as a function of visual magnitude for the three clusters 
taken together, over the range —4.1 < V-V+ om < —2.0 mag, |5(V-J)| < 0.07 mag, where 
S( V-I) is the horizontal distance of a given point from the best-fitting parabola. This 
region is denoted by the curve-sided box in Fig. 17. Then the net offset of the data for each 
Table 4 individual cluster from the fitted parabola was determined, as listed in Table 4: the table 
gives (column 2) the mean differential offset for the giant branch of each cluster, (column 3) 
the median offset among the giants for each cluster, (column 4) the standard deviation 
about the mean horizontal offset for the each cluster, (column 5) the number of stars 
contained in the box, and (column 6) the median value of a( V-I) among the stars in the 
box. We note that the perceived standard deviation is quite a bit larger than our estimates 
of the photometric errors for the individual stars. There can be several reasons for this: 
(a) we may have underestimated our photometric errors; (b) there could be actual scatter 
among the giants in each cluster; (c) the giant branch might not be a perfect parabola over 
this range of magnitude; and (d) the derived standard deviation could be dominated by the 
random errors of stars with individual er( V-I) values larger than the median value. 

Using the isochrones described above, we found that at [Fe/H]^ —1.41 and age 
~ 12 Gyr, the relationship between RGB color offset and age is S[A(V-I)]/St = —0.014 
mag/Gyr. If we take the relative giant-branch offset of Pal 3 and Pal 4 to be of order 
0.000 ± 0.003 mag, then we would argue that these clusters are coeval to a level of about 
0.2 Gyr if the difference in their abundances can be neglected (although we will argue below 
that the abundance differences should not be neglected). On the other hand, the Eridanus 
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giant branch appears to be offset from the giant branches of the two Palomar clusters by 
some 0.01 ± 0.004 mag, which would suggest that Eridanus is younger than the other two 
clusters by some 0.7 ± 0.3 Gyr under the same set of assumptions. 

It is also quite apparent in Fig. 17 that by aligning the three cluster sequences at 
the main-sequence fiducial magnitude V+ om , we have also brought the three clusters' HBs 
into excellent coincidence. The small number of HB stars in each cluster perhaps makes it 
dangerous to try to be overly quantitative (note in particular that the two double circles 
on the HB bluer than ( V-I) — { V-I)to — 0-12 represent mean values for two of the Pal 3 
RR Lyrae variables; the other five were excluded by the a {V-I) < 0.05 criterion), we feel 
that we may conservatively estimate that the AV^p values for the three clusters do not 
differ by more than 0.05 mag or so (we present actual measured values in Table 6 below). 
As our isochrones suggest that at [Fe/H]?s —1.41 and age 12 Gyr the turnoff fades by 
some +0.085 mag/Gyr in V, this would suggest an upper limit of some 0.6 Gyr for the total 
age spread among these clusters if their abundances are similar. On the other hand, it is 
commonly believed that at a fixed age the V^-band luminosity of the HB decreases by an 
amount of order 0.2 mag for every 1 dex increase in the overall heavy-element abundance 
{e.g., Dorman 1993; Chaboyer, Demarque & Sarajedini 19966; Gratton et al. 1997; our 
theoretical models are consistent with this: VandenBerg et al, in preparation). If this is 
true, and if Pal 3 is of order 0.2 dex more metal poor than Pal 4 and Eridanus, then its HB 
should be of order 0.04 mag more luminous than in those other clusters. As we will discuss 
below, we believe that an offset of this order of magnitude is within the measuring error of 
the HB magnitude. 
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3.2. Comparison with the Near-Halo Clusters M5 and M3 

With the CMD data for these three clusters we are in a position to begin answering a 
two-decades-old question: are the outer halo clusters with anomalously red HB morphologies 
truly younger than their counterparts in the nearby halo? Even with these data from 
HST, the answer to this question will depend on the estimates of the abundances of the 
outer halo clusters as well as for our representatives of the nearby halo, M3 and M5. (M5 
is currently inside the 8kpc boundary that nominally separates the inner halo from the 
intermediate halo, but its large proper motion — Cudworth & Hanson 1993 — indicates 
that this situation is only temporary.) From the arguments presented in the Appendix, 
we believe that it is adequate at present to consider that Pal 4, Eridanus, and M 5 have 
essentially the same metallicity, while the metal abundances of Pal 3 and M 3 are similar to 
each other and perceptibly lower than those of the other three clusters. 

While the principal sequences of the intermediate-halo clusters M 3 and M 5 can be 
compared to those of the outer-halo clusters by the same techniques as we have just applied 
to the comparison of the outer-halo clusters among themselves, it is best to remember that 
there are some underlying differences. The comparisons among the outer-halo clusters are 
direct: the observations were made with the same telescope (HST) and camera (WFPC2), 
and were subjected to an identical analysis and calibration procedure. M 3 and M 5, 
in contrast, were each observed with several different ground-based telescopes and on 
numerous occasions. Still, in many cases both clusters were observed on the same nights, 
and all observations have been referred to a common photometric system with the utmost 
rigor currently possible. However unlikely, there remains the possibility of some undetected 
difference between the HST and ground-based photometric systems that will render a 
comparison of the nearby halo clusters to the outer-halo ones less robust than comparisons 
that are strictly internal to either data set. Uncertainties in the zero points of the color 
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and magnitude scales are not of immediate concern; what might cause a significant problem 
would be uncorrected color dependences in the transformations from the instrumental to the 
standard magnitude systems, or a nonlinearity in the magnitude scale of either the space- 
or ground-based photometry. We stress that at present we have no evidence that such 
problems exist. Nevertheless, it would be wise to retain an extra quantum of skepticism in 
the backs of our minds as we interpret the comparison of the HST data to the ground-based 
observations. 

We have repeated the determination of the relative giant-branch shifts for the five 
clusters just as in the previous section. This time we combined the rectified (V, V-I) data 
for all five clusters and refit the mean giant-branch parabola using the same constraints as 
before: —4.1 < V^-V+o.os < —2.0, | 8{V-I)\ < 0.07. We were a little more stringent in our 
selection of stars for M3 and M5, retaining only those stars with individual color errors 
a {V-I) < 0.03 mag. After the parabola had been fit, we determined the net offsets and 
scatter for each of the individual clusters relative to the mean giant branch of all of them 
Table 5 as before; the results are listed in Table 5. The implications are that, if Pal 4, Eridanus, 
and M5 all have the same composition — and we stress again that this means both [Fe/H] 
and [cc/Fe] - - then Pal 4 and Eridanus are younger than M5 by some 1.4 and 2.1 Gyr, 
respectively. The difference between Pal 3 and M 3 — again, assuming that these two 
clusters have the same chemical abundances — is similar: the net relative giant branch 
offset would imply that Pal 3 is younger than M3 by some 1.7 Gyr. 

The alternative technique of determining the relative ages of clusters from the 
magnitude differences between their HBs and their main-sequence turnoffs requires a 
robust, consistent measurement of their HB magnitudes as well as of their turnoffs. The 
upper- main-sequence and subgiant-branch fiducial points given in Table 3 serve the latter 
purpose; now it remains to try to estimate the fiducial HB magnitudes with similar rigor. 



-22 - 



We have computed zero-age horizontal branches (ZAHBs) consistent with the predictions of 
these main-sequence and giant-branch evolutionary tracks. We have employed them much 
as we used the isochrones in §3, simply as drafting splines of approximately the right shape 
to match the shape of an actual cluster HB. For comparison with Pal 3 and M 3, we have 
adopted theoretical ZAHB sequences with [Fe/H] = -1.54, [a/Fe] = +0.3, and Y = 0.2362; 
Pal 4, Eridanus, and M5 were compared to theoretical curves with [Fe/H] = -1.41, [a/Fe] 
= +0.3, and Y = 0.2366. Unlike the case with the turnoff isochrones, in this case we have 
allowed no freedom for horizontal translation of the theoretical curves, but have adopted 
their predicted colors and estimated reddening values verbatim, and have performed vertical 
shifts only in matching the theoretical ZAHBs to the lower envelope of the dereddened 

Fig. is cluster data. Figure 18 shows these comparisons. The apparent visual magnitude of the 
fitted ZAHB was read out at an intrinsic color of (V-I) = 0.60, and the resulting Vhb 

Table 6 values are listed in Table 6. The external accuracy of these derived HB magnitudes is 
clearly dependent upon the requirement that the theoretical models accurately reproduce 
the luminosity variation of the HB with color, both at the red end (for the three outer-halo 
clusters), and at the blue end (for the two nearby clusters): if these curves are seriously 
wrong in detail, errors in the vertical placement of the ZAHB of several hundredths of a 
magnitude could easily result. 

The precison of the HB magnitudes is also plainly jeopardized by the small number 
of HB stars in each cluster. For the three distant systems, there is little more that can 
be done; these diagrams contain most of the HB stars that Nature has provided. For 
M5, at least, we have a convenient check. Sandquist et al. (1996) have published (V,I) 
photometry for a large sample of stars in this cluster. Johnson & Bolte (1998) find some 
modest errors in the zero points and possibly the scale of the Sandquist et al. magnitudes 
and, based on our own analysis, we agree. However, considering only those stars common 
to the Sandquist et al. sample and our current data set with 14.5 < V < 16.5 — 60 
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stars in all — we find net differences of AV = +0.009 ± 0.0029 mag (standard error of 
the mean, ±0.022 standard deviation of one difference) and Al = —0.006 ± 0.0026 mag 
(s.e.m., ±0.020, s.d.) in the sense (present - Sandquist et al). The data of Sandquist et 
al. for stars with cr(V-I) < 0.03 mag and lying more than 96" from the center of M5 are 
Fig. 19 plotted, after correction for these offsets, in Figure 19 along with the same fitted ZAHB as 
before. [Note that here (m-M)y does not necessarily mean —5 ± 51ogrf(pc) ± Ay', rather it 
means my (observed)— My (predicted).] It seems that in this case, at least, our derived HB 
magnitude is consistent with the larger data set. 

The last column of Table 6 lists the values of Vpo — Vhb that result from the difference 
between these HB magnitudes and the turnoff magnitudes derived from the isochrone fits. 
We see that for the three outer-halo clusters, the magnitude differences between the HB 
and the turnoff are nearly indistinguishable. Likewise, the difference between M 3 and M 5 
is infinitesimal, but there is a difference of approximately 0.15 mag in the turnoff-to-HB 
magnitude difference (or the subgiant-branch-to-HB difference) between the nearby and 
remote clusters. In the approximation where all clusters have effectively the same chemical 
abundances, this would imply that the three outer-halo clusters are all roughly 1.5 - 2 Gyr 
younger than the two near-halo clusters. 

To further quantify this conclusion, it is necessary to look more closely at the 
predictions of the current set of isochrones, in particular at the dependence of A(V-I) 
and AV" on chemical abundance in comparison to their dependence on age. The solid 
Fig. 20 curve in Figure 20 shows the predictions of VandenBerg's isochrones for the evolution 
of the turnoff-to-giant-branch color difference (horizontal axis) and the turnoff-to-HB 
magnitude difference (vertical axis) as functions of age, from 8 to 18 Gyr, for abundances 
of [Fe/H] = -1.41 and [a/Fe] = ±0.3. The horizontal tick marks intersect the curve at 
the model predictions for incremental age increases of 2 Gyr. The arrow illustrates the 
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effect on the models produced by an assumed increase of +0.2 dex in either [Fe/H] or 
[O/Fe], according to this particular set of isochrones. There is considerable uncertainty 
in the horizontal placement of the solid curve in this diagram, as the predicted colors 
are dependent both upon uncertain details of convection theory and upon the assumed 
conversion from theoretical effective temperatures to observational colors. The slope of 
the curve could also be somewhat in error if the assumed dependence of convection or the 
temperature-to-color conversion on temperature itself is in error in a seriously non-linear 
fashion. Similarly, there could be some uncertainty in the absolute vertical positions of the 
tickmarks if the predicted HB luminosities are in error, or if the conversion of bolometric 
magnitudes to V-band magnitudes is incorrect for either the turnoff or HB stars. Therefore 
we particularly wish to de-emphasize any implication that we are measuring absolute ages 
for any of these clusters, at least until some of the remaining uncertainties in the model 
physics have been resolved. However, the size of the intervals between the vertical positions 
of the tick marks, representing increments of order 15-20% in age, should be relatively 
secure. The locations in this plane of the three outer-halo clusters are represented by open 
symbols, a triangle for Pal 3, a square for Pal 4, and a circle for Eridanus, while the two 
nearby clusters are indicated by filled symbols: a triangle for M 3 and a pentagon for M 5. 
Representative one-sigma error bars of ±0.005 mag in A(V-I) and ±0.05 mag in AV - 
due primarily to the uncertain placment of the HB lower envelope — are illustrated in the 
bottom right corner. 

A consistent interpretation of the two age indicators can be made for the five clusters. 
// we assume that the theoretical isochrones employed here correctly predict the dependence 
of AV and A(V-I) on age and chemical abundance, then: 

1. Those same models slightly (i.e., by about 0.017mag) overestimate the color difference 
between the giant branch and the turnoff, and the solid curve should be shifted as a 
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body to the left. 

2. The data are consistent with an age difference of some 2 Gyr between M3 and Pal 3 
(Pal 3 younger) if their abundances are the same. Similarly, Pal 4 is of order 1.5 Gyr 
younger than M 5 if their abundances are identical. 

3. The difference between the CMDs of Pal 4 and Eridanus could be due to 

a. an age difference of somewhat less than 1 Gyr (Eridanus younger) if their 

abundances are the same; or 

b. an abundance difference of some 0.2 dex with Eridanus being the more metal poor, 

in which case the age difference would be small. 

The available spectroscopic data suggest that Eridanus may be some 0.14 dex more 
metal-poor than Pal 4, so the evidence for an age difference between the two is 
tenuous. 

4. Even though the point for Pal 3 coincides with that for Pal 4 in this diagram, it is 
likely that Pal 3 is actually ~ 1 Gyr older than Pal 4, because its point will have 
been displaced upward and to the right in this diagram as a result of its lower metal 
abundance. 

5. Similarly, M5 is almost certainly more metal-rich than M3, so by the same argument 
as in point 4 above, M5 must also be younger than M3, probably by of order 1 Gyr 
or somewhat over. 

4. Discussion 

The two relative age indicators, A( V-I) and AV, suggest a single self-consistent 
interpretation, namely that all three outer halo clusters are younger than their near-halo 
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counterparts by 1.5 to 2 Gyr, if we have correctly estimated the clusters' chemical 
abundances. Sarajedini's (1997) result for Pal 14 places this fourth outer-halo cluster in the 
same age ballpark (with, of course, the same caveats regarding uncertainties in the chemical 
abundance ratios and poorly understood processes in stellar interiors). Based on LDZ's HB 
star models, if the red HB of Pal 3 is due to an age difference between this cluster and the 
average [Fe/H]~ —1.5 cluster in the near halo, Pal 3 is ~2 Gyr younger than its near-halo 
counterparts. For Pal 4 and Eridanus, with HB types (B — R)/(B + V + R) = — 1, Lee's 
method can assign only a lower limit of > 2 Gyr to the age difference, since the models 
for all ages younger than this imply wholly red HBs. These differential ages based on HB 
morphologies are in the same sense and of approximately the same size as those implied by 
our turnoff and subgiant comparisons. 

Another angle from which to view this result is to ask, "What other parameters of 
the outer-halo clusters would have to be different if they are to have the same ages as 
the near-halo clusters?" We conclude that their age differences would be reduced to less 
than ~ 1 Gyr (below which we would find it difficult to argue convincingly for genuine 
age differences) if we have overestimated the heavy-element abundances for the outer-halo 
objects — either [Fe/H] or [cc/Fe], or a combination of the two — by 0.3 dex or more. We 
cannot rule out this possibility or, for that matter, less frequently discussed ones such as 
differences in helium abundance, mass loss rates on the giant branch, or stellar rotation. 
As we discuss in the Appendix, the [Fe/H] determinations for the outer-halo clusters 
are considerably weaker than those of the near-halo ones, and we have no quantitative 
information at all on their [a/Fe] values. It is intriguing in this context that recent 
spectroscopic observations by Brown, Wallerstein & Zucker (1997) indicate that the young 
halo globular clusters Ruprecht 106 and Palomar 12 have relative oxygen and general 
a-element abundances [O/Fe] and [ct/Fe] ~ 0.0 to +0.1, in contrast to the value +0.3 
to +0.4 observed in nearby, "normal" halo objects, although in the case of oxygen it is 
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unknown whether the abundance ratio in Rup 106 and Pal 12 is primordial or has been 
altered by mixing. If Pal 3, Pal 4, and Eridanus were to have primordial abundances [a/Fe] 
~ 0.0, then by basing our metallicity on measurements of Ca calibrated via nearby clusters 
with [a/Fe]> we have assigned the wrong [Fe/H]. The sense of this error is such that 
it would reduce the inferred age difference between them and the near-halo clusters. We 
therefore believe it is of the utmost importance to obtain high-resolution spectroscopic data 
for the outer-halo clusters. 

We have been scheduled observing time on HST to obtain color-magnitude diagrams 
for the last two remaining globular clusters beyond 50kpc from the Galactic center, 
Palomar 14 and AM I. At the same time, other research teams are using HST to pursue the 
study of the globular clusters near the Galactic center, and the dwarf galaxies that share the 
outer halo with the Palomar-like clusters, while still others are improving the available data 
for the near-halo objects. Thus, this is probably not the time to draw definitive conclusions 
regarding the early evolution of the Galactic spheroid. We can nevertheless draw some 
provisional conclusions that may be tested and rejected or refined by future work. As 
always, these conclusions are subject to the provisos that we have correctly estimated the 
chemical abundances of the objects under discussion and that no important physics has 
been omitted from the differential predictions of stellar evolution; simple errors in distance 
scale, convection theory, opacities, elementary particle physics, and anything else can be 
neglected provided they affect all globular clusters comparably, and provided we restrict 
ourselves to investigating the order and relative time intervals over which things happened, 
as opposed to the exact moment that any one event took place. 

These results make a good case for age as a second parameter: the mean age of the 
outer-halo clusters with anomalously red HBs appears to be lower by ~ 1 — 2 Gyr than 
that of the nearby, "normal" clusters. Although the original motivation for many programs 
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investigating ages for GGC has been to decide the mix of ELS and SZ that is required to 
match the cluster age distribution, it is clear that the GGC age distribution alone does not 
lead directly to a definitive formation scenario for the Galaxy. A range in GGC ages (> a 
few Gyr) is clearly inconsistent with the most naive version of a monolithic collapse model, 
but a general large-scale lumpy collapse accompanied or followed by various degrees of infall 
or accretion are allowed. Hierarchical formation models, like SZ, can more easily accomodate 
a range of ages than ELS, but they do not make specific predictions of the amount of time 
required to assemble the Galactic Halo and they do not predict an age distribution for 
GGC. Observationally, the Sagittarius dwarf spheroidal galaxy provides a striking example 
of a "late" accretion event that will add field stars and globular clusters to the Galaxy, yet 
not significantly alter the age dispersion of the GGC system. Still, it would be premature 
to conclude on the basis of the Sagittarius example that accretion was necessarily the sole, 
or even the dominant, mode of halo formation. Based on high-resolution simulations of the 
formation of galaxies after recombination (e.g., Sommerville & Primack 1998), there is a 
large overlap in the progenitor collapse times of objects that become galaxies of different 
masses. Objects that eventually become LMC-like galaxies can first start to form stars 
before, at the same time or after Milky Way-like objects. The age distribution is a crucial 
piece of the picture, but the final story will depend on other information such as detailed 
abundance patterns and kinematics for GGC halo populations (e.g., Carney et al. 1997; 
Zinn 1993). 

In contrast to theories that attempt to describe the hierarchical formation process 
of the Galaxy as a whole, the Harris & Pudritz (1994) formation model for the globular 
clusters themselves does make specific quantitative predictions for the amount of time 
required to assemble these systems within the outer Galactic halo. Far out in the newly 
forming Galaxy, the gas clouds (supergiant molecular clouds or SGMCs) in which the 
protoclusters are presumed to have been built would have had characteristically larger 
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linear sizes, to match the lower ambient pressure that was found there, in comparison to 
the denser conditions in the rapidly collapsing Galactic core. The timespan over which 
protoclusters can form scales directly with the internal free-fall time of the SGMC; thus, the 
cluster formation timespread will increase as the surrounding pressure drops. If the SGMC 
lifetime in the mid-halo region was typically ~ 0.5 Gyr (see Harris & Pudritz equation 5.9 
and related discussion), then it could have been > 2 Gyr in the outermost regions. 

Finally, as already noted by van den Bergh (1998), the GGC which have been shown 
to be young are all at the low-luminosity end of the My distribution and are all beyond 
Rqc = 15 kpc. The massive clusters beyond this Galactocentric radius do appear to be 
co-eval with the oldest near-halo clusters. Van den Bergh interprets this to indicate that 
the mass of newly-forming clusters decreased with time in the outer halo. Our present data 
agree with this interpretation. It will be very interesting to determine reliable ages for the 
remaining clusters at large Rqc (all relatively low-L objects) and to pursue similar studies 
of the low-L GGC closer to the center of the Galaxy. 
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A. [Fe/H] and [a/H] Determinations from the Literature 

Extracting reliable relative or absolute age information from CMDs depends 
upon having accurate relative and absolute values for the chemical compositions of 
their constituent stars. In particular, the morphology of the HB is sensitive to (often 
unmeasurably) small differences in many parameters, among them age, [Fe/H], [CNO/Fe], 
[a/Fe], Y, mass, core rotation, and so on. Thus, constraining the role played by age 
differences in producing a range of HB morphologies among clusters requires understanding 
the roles played by all other variables. A special emphasis emerges on understanding 
possible chemical differences which — unlike, say, interior rotation — we can, in principle, 
directly measure. (However, note recent papers by Rutledge, Hesser & Stetson 1997 and 
King et al. 1998, which demonstrate that there remains considerable doubt about the 
significance of the chemical abundances that can be inferred from existing spectroscopic 
observations.) Similarly, differential comparisons among remote and nearby halo objects 
offer a powerful tool to constrain possible age gradients as a function of Galactocentric 
radius, provided that we can identify near-halo clusters with compositions comparable to 
those of the outer-halo objects. In this section, we review information from the literature 
on our outer halo clusters and the comparison clusters we have chosen for them. 

A.l. General remarks 

It is a risky undertaking to combine abundances from different literature studies 
because there is not a general agreement on the overall abundance scale. The most widely 
used general metallicity scale is that of Zinn & West (1984; hereinafter ZW84), which 
is ultimately based upon photometric, spectrophotometric, and photographic-spectrum 
abundance determinations from the 1970's. Although some recent studies have called into 
question the zero point and scaling of absolute metallicity values from ZW84 {e.g., Carretta 
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and Gratton 1997), it appears to provide generally reliable relative rankings of cluster 
metallicity (e.g., Rutledge et al. 1997). Of greatest relevance to our present study, we note 
that in the range —1.0 < [Fe/H] < —1.9, the Zinn-West values are systematically lower 
by 0.23 dex than the Carretta-Gratton values. Similar trends are seen when comparing 
Rutledge et aVs large, homogeneous survey of Ca II triplet measures to the two scales. We 
are fortunate that Pal 3, Pal 4 and Eridanus have all been the subject of abundance studies 
which allow them to be differentially compared to our two "standard" clusters M 3 and M 5 
(see details below); to some extent, it is this relative ranking by metallicity that is most 
important for determining the range of ages among the clusters. 

Because we will put high weight on [Fe/H] determinations for the outer-halo clusters 
that have been inferred from the lines of Ca in the near-IR, we need to be concerned not 
only with [Fe/H] but also with the size of the cluster-to-cluster variations in [Ca/Fe] or more 
generally [a/Fe]. In recent years considerable observational evidence has accumulated that 
near-halo clusters and field stars generally display roughly constant ratios of a to iron-peak 
elements, with the former being more abundant than the latter by several tenths of a dex 
for [Fe/H]< —1.0; the a-to-iron ratio approaches the solar value as [Fe/H] increases from 
-1.0 to 0.0, at least among the field stars surveyed. In a review of elemental globular cluster 
abundances based again on modern high- dispersion studies, Carney (1996) concludes that 
there is no variation of [a/Fe] as a function of metallicity for 14 nearby clusters, all with 
[Fe/H] < -0.6 to -0.7. He further summarizes the evidence for an a-element enhancement 
relative to iron of [a/Fe] ~ +0.28 dex for these objects. 

The great distances of the most remote clusters make fine abundance analyses 
challenging, and there is so far little direct evidence bearing on [a/Fe] trends in the 
outer halo. Indeed, there are only a few such studies for clusters at intermediate 
Rgc- Brown, Wallerstein & Zucker (1997) have found that two giant stars in each of 



-32 - 



Ruprecht 106 and Palomar 12 exhibit solar ratios of the a elements (specifically, Mg, Si, 
Ca, and Ti), yet the clusters' overall heavy-element abundances ([Fe/H]=— 1.45 ± 0.10 and 
— 1.0 ± 0.10, respectively) fall in the regime where most halo objects exhibit enhanced a 
elements. Moreover, for Rup 106 they detected an oxygen anomaly, in that its value of 
[0/Fe]= 0.0 ± 0.13 is about 3a below that normally found for stars with [Fe/H] ~ —1.5. 
Could such anomalies also occur in other outer-halo clusters? We have to be cautious in 
interpreting the latter measurement because [O/Fe] has been shown in several clusters 
to decrease with increasing RGB luminosity, suggesting that CNO-processed material is 
being dredged into the star's atmospheres. The low [O/Fe] for the brightest Rup 106 stars 
may therefore not reflect the initial oxygen abundance of the stars. Based on six giants in 
NGC 7006 which, at R GC = 38kpc, is half of the way out to the far halo, Kraft et al. (1998) 
find [Ca/Fe]~ +0.2 — consistent with the general trend of the near-halo clusters but also, 
at the other extreme of the error bar, possibly consistent with less a-element enhancement 
than usual in the near halo. With its very large space velocity (Cudworth & Hanson 1993), 
M5 might itself be considered an outer halo denizen that shows "normal" [a/Fe]. 

For Population II field stars, there have been additional indications that [a/Fe] may 
be solar or subsolar for some objects with very large space velocities and therefore large 
inferred apogalactic distances (e.g., Carney et al. 1997; King 1997). It is too soon to draw 
any general conclusions about systematic trends in [a/Fe] with Rqc- m what follows we 
must keep in mind that these remote objects may have [a/Fe] ratios that are less enhanced 
than those of our fiducial clusters M3 and M5, and from the other calibrating clusters for 
the metallicity scales. 
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A.2. Eridanus 

Of the three outer-halo clusters, only Eridanus was included in the original ZW84 
compilation; they assigned it [Fe/H]= —1.4, the same as M 5 and ~ 0.25 dex more metal-rich 
than M3. Armandroff & Da Costa (1991) obtained Ca II triplet spectra for two Eridanus 
giants along with numerous giants in six nearby well-studied clusters. Their reduced calcium 
equivalent widths suggested that Eridanus falls about one-third of the way from NGC 6752 
([Fe/H]= -1.55) to NGC 1851 ([Fe/H]= -1.26), two of their "standard" clusters. Thus 
from the Ca II measures alone, they derived —1.41 ± 0.11 on the ZW84 scale. (The same 
placement relative to NGC 6752 and NGC 1851, interpreted through the Rutledge et al. 
version of ZW84 and the Carretta-Gratton scales would be —1.44 and —1.14, respectively. 
All three of these estimates are internally consistent to within ±0.1 dex, given the 0.23 dex 
offset between ZW84 and CG97 mentioned above.) Armandroff & Da Costa combined their 
Ca II abundance values with two earlier estimates (—1.50 ±0.15 by Da Costa & Armandroff 
(1990), derived from the position of the Eridanus giant branch in the (Mj, (V— I) Q ) plane; 
and Da Costa's (1985) value of [Fe/H]= -1.35 ± 0.2 determined from (B-V) , s ) to adopt a 
mean —1.42 ± 0.08 dex for the cluster (again on the ZW84 scale). 

Earlier, Ortolani & Gratton (1989) had analyzed low-signal-to-noise, low-dispersion 
spectra of two Eridanus giants and four Pal 3 giants (see below). They concluded 
that Eridanus falls between NGC 3201 and NGC 4590 in metallicity and assigned it 
[Fe/H]= —1.6 ± 0.3 on their abundance scale; since ZW84 assigns metallicities of -1.6 and 
-2.1 to NGC 3201 and NGC 4590, respectively, Ortolani & Gratton's spectra might imply 
[Fe/H] ~-1.8 on the ZW84 scale for Eridanus, but the error bar is large. 
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A.3. Pal 3 

Pal 3 and Pal 4 were included in the investigation of Armandroff, Da Costa & Zinn 
(1992) in which abundances were measured from the Ca II triplet. Spectra for three Pal 3 
giants were measured and the reduced equivalent widths placed Pal 3 more metal-poor than 
M5 and similar to M 13. In this study, Pal 3 was also ranked more metal-poor than Pal 4. 
On the ZW84 scale, Pal 3 was assigned [Fe/H]= -1.57 ± 0.19. 

Pal 3 was also part of the Ortolani & Gratton (1989) study mentioned above. From 
low-dispersion spectra of four giants, Pal 3 was assigned [Fe/H]= —1.6 ± 0.3 (the same as 
Eridanus from that study) although the star-to-star scatter was large. From their CMD for 
Pal 3, Ortolani & Gratton estimated E(B-V) = 0.03 ± 0.02 and (B-V) 0tg = 0.69 ± 0.04, 
leading to [Fe/H]= — 1.88±0.15. They adopted — 1.7±0.2 from this combination of spectral 
and CMD indicators. 

A. 4. Pal 4 

Ca II triplet measures for three giants led Armandroff et al. (1992) to conclude that 
[Fe/H]= —1.28 ± 0.20 for Pal 4 on the ZW84 scale, compared to earlier estimates from 

CMDs of [Fe/H]= -1.7 (Christian & Heasley 1986) and [Fe/H] 1 (Reed & Harris 1986). 

They further note that their triplet metallicity determination would be —1.35 ± 0.14 were 
NGC 6171 excluded from their calibration. ^From their Figure 4, Pal 4 appears to be clearly 
more metal-rich than M 13 and Pal 3 and comparable to or slightly more metal-rich than 
M5. 
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A. 5. Halo Cluster Metallicity Summary 

The previous studies are in good general agreement in finding the Eridanus metallicity 
to be ~ -1.4 on the ZW84 scale. To the extent that NGC6752 is similar in [Fe/H] to 
M3 and M 13, Eridanus appears to be more metal-rich than M3. Based primarily on the 
Armandroff et al. Ca II triplet work, we conclude that Pal 3 is closer to M 3 and M 13 
in metallicity than it is to M5. Furthermore, by appealing to the differential studies of 
Armandroff et al. (1992) and Ortolani & Gratton (1989) we suggest that Pal 3 is the most 
metal-poor of our three clusters. 

For purposes of comparison, our adopted [Fe/H] abundances for all five clusters, on 
the ZW84 scale, are listed in Table 6. The values for M3 and M5 are taken from the 
compilation of Harris (1996); each of the five is likely to be externally uncertain at the 
±0.15 dex level at least (see Paper I). Nevetheless, they indicate formally that there is no 
significant difference between M3 and Pal 3, or between M5 and Pal 4. On the strength 
of both the metallicity measurements and the CMD comparisons shown above, we believe 
Eridanus to be most nearly similar to M 5 as well. 
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Fig. 1. — Color- magnitude diagram for the outer-halo globular cluster Pal 3, derived 
from our WFPC2 photometry with the HST. Only stars with photometric standard errors 
<r( V-I) < 0.10 mag have been plotted, and data from all four CCD fields have been combined. 
The circled cross at the blue end of the HB represents the mean magnitude and color of 7 
RR Lyrae stars in the cluster. 

Fig. 2. — Color- magnitude diagram for Pal 4, as in Fig. I. 
Fig. 3. — Color- magnitude diagram for Eridanus, as in Fig. 1. 

Fig. 4. — Ground-based (V, V-I) color-magnitude diagram for the nearby globular cluster 
M3. Only stars with photometric standard errors cr(V-I) < 0.10 mag have been plotted. 

Fig. 5. — Ground-based (V, V-I) color-magnitude diagram for the nearby globular cluster 
M5, as in Fig. 4. 

Fig. 6. — Light-curve fragments for seven candidate RR Lyrae variables in the outer- 
halo globular cluster Pal 3, based upon our HST observations. Open squares represent 
instantaneous /-band magnitudes, and closed circles represent V-band magnitudes. 

Fig. 7. — The open squares are our fiducial points for the CMD of Pal 3, as determined by 
mean colors after three iterations of cr-clipping in bins of height 0.2 mag in V. The solid curve 
is a robust maximum-likelihood fit of which isochrone that best matches the observed CMD 
from among a wide variety of isochrones that were tried. For the HB level and brighter, the 
individual HB and RGB stars are plotted (crosses). 

Fig. 8. — Fiducial points and isochrone fit for the CMD of Pal 4, with symbols as described 
in Figure 7. 

Fig. 9. — Fiducial points and isochrone fit for the CMD of Eridanus, with symbols as 
described in Figure 7. 
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Fig. 10. — Fiducial points and isochrone fit for the CMD of M3, with symbols as described 
in Figure 7. 

Fig. 11. — Fiducial points and isochrone fit for the CMD of M5, with symbols as described 
in Figure 7. 

Fig. 12. — Theoretical isochrones for chemical-abundance ratios [Fe/H] = -1.41 and [a/Fe] 
= +0.3, and ages of 8 to 16 Gyr in steps of 2 Gyr, illustrating the predicted dependence of 
CMD morphology on age, as described in the text. 

Fig. 13. — Theoretical isochrones for chemical-abundance ratios [Fe/H] = -1.14, -1.41, and 
-1.61; [a/Fe] = +0.3; and an age of 12 Gyr, illustrating the predicted dependence of CMD 
morphology on relative iron abundance. The isochrones have been shifted horizontally and 
vertically as in Fig. 12. 

Fig. 14. — Theoretical isochrones for [Fe/H] = -1.41; [a/Fe] = 0.0, +0.3, and +0.6; and 
an age of 12 Gyr, illustrating the predicted dependence of CMD morphology on relative 
a-element abundance. The isochrones have been shifted horizontally and vertically as in 
Fig. 12. 

Fig. 15. — Fiducial points for the three outer-halo clusters (Pal 3, Pal 4, and Eridanus), 
dereddened according to the values listed in Table 6, and shifted vertically to match their 
red HBs (the short horizontal line at My = 0.70 represents the RHB level of each cluster). 
The distance scale for My(HB) has been chosen arbitrarily for comparison purposes. Note 
that the mean points for Pal 4 and Eridanus are essentially indistinguishable; the giant 
branch of Pal 3 lies ~ 0.02 magnitude to the blue of the GBs of the other two clusters. 

Fig. 16. — The mean points for Pal 3, Pal 4, and Eridanus have been re-registered as 
prescribed by VandenBerg, Bolte, & Stetson (1990): the CMDs have been shifted horizontally 



-43- 



so that the MSTO colors coincide (vertical line), and vertically so that the point on the 
main sequence 0.05 mag redder than the MSTO coincides (horizontal line). The horizontal 
shifts were taken from the isochrone fits (Table 3), which produce estimated turnoff colors 
systematically 0.007 mag bluer than a-clipped normal points given in Table 2; thus there 
appears to be a slight offset of the turnoff points from A( V-I) = 0.00. This offset is nearly 
identical for the three clusters, so no material harm is done to the differential analysis. For 
a given composition, an age difference appears as an offset in the RGB color. Pal 4 and 
Eridanus appear to be indistinguishable in age, while Pal 3 may be slightly older. If the 
clusters do not have identical abundances, the interpretation of this diagram becomes more 
complex, as described in the text. 

Fig. 17. — The data for Pal 3, Pal 4, and Eridanus have been registered as in Fig. 16, 
except in this case the stars have been plotted individually. The curve-sided box represents 
the region where the relative displacements of the three giant branches were determined, as 
listed in Table 4. Only stars with a {V-I) < 0.05 mag have been plotted. Symbol types are 
as in Fig. 16. 

Fig. 18. — Stars in the region of the red HB have been plotted for (from top to bottom) 
Pal 3, Pal 4, Eridanus, M3, and M5. The photometry has been dereddened according 
to the color excesses in Table 6, with the assumption that E{V-I) = 1.3E{B-V). The 
solid curves represent our eyeball fits of theoretical ZAHB loci to the lower envelope of the 
horizontal-branch stars for each cluster. 

Fig. 19. — Horizontal-branch region of the CMD of M 5 based upon the V, I data of Sandquist 
et al. (1996), after correcting for minor zero-point shifts in the photometry as described in 
the text. The solid curve is the same ZAHB locus as in the bottom panel of Fig. 18. Here, 
{m-M) v denotes the vertical shift required to register theoretical luminosities M v with 
observed V-band magnitudes; due to uncertainties in the modelling, this may not accurately 
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represent the actual geometric distance of the cluster. 

Fig. 20. — The relationship between the giant-branch-to-turnoff color difference (abscissa) 
and the horizontal-branch-to-turnoff magnitude difference (ordinate), according to the 
present set of isochrones. The solid curve represents the predicted track followed by a 
globular cluster with [Fe/H] = -1.41 and [a/Fe] = +0.3 as it ages from 8 Gyr (top right) 
to 18 Gyr (bottom left) according to the present set of theoretical models. The absolute 
placement of the curve and tick marks are uncertain due to incomplete knowledge of the 
physics of stellar interiors, but the size of the intervals, representing age differences of 2 Gyr, 
or 15-20%, should be correct in a relative sense. The arrow represents the amount by which 
the models would be shifted for an iron or a-element abundance 0.2 dex larger. The positions 
in this diagram currently occupied by the outer-halo globular clustars are represented by 
open symbols: triangle (Pal 3), square (Pal 4), and circle (Eridanus). The near-halo clusters 
are represented by closed symbols: triangle (M3) and pentagon (M5). The cross in the 
lower right represents our estimate of the random observational errors (±lcr) in the two 
coordinates. 
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Table 1. Notional Periods for Candidate RR Lyrae Variables 
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Table 2. Normal Points for Program Clusters 
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Table 2 — Continued 
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n qi 7 _L n ooft 
U.ol ( zt u.uuo 


99 
ZZ 


9n 

ZU 


1 Q 

iy 


93 ^ns -i- n nn.fi 
zo.ouo m u.uuo 


n 70 A J- OOfi 
U. ( Z'i ZIZ u.uuu 


R7 
O / 


uo 


uo 


23.714±0.007 


0.644 ±0.002 


61 


58 


55 


23.901 ±0.007 


0.617 ±0.003 


70 


67 


66 


24.098 ±0.006 


0.621 ±0.002 


96 


93 


89 


24.299 ±0.006 


0.620 ±0.003 


109 


103 


99 


24.499 ±0.005 


0.632 ±0.003 


132 


129 


124 


24.699 ±0.005 


0.645 ±0.002 


159 


151 


143 


24.908 ±0.004 


0.665 ±0.003 


167 


160 


153 


25.095 ±0.005 


0.682 ±0.003 


165 


161 


151 
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V 


V-I 


ni 


^2 


n 3 


9^ "?nn -i- n nnA 

ZO.OUU in U.UU4: 


n 7nfi 4- n nn9 

U. ( Uu It U.UUZ 


00 A 
ZZ^ 


90Q 

zuy 


iy 


9^ ^nn -i- n nn/i 

ZO.OUU it U.UU4 


n 79R -i- n nn*? 

U. 1 ZO it U.UUO 


99Q 

zzy 


91 ^ 
Z10 


OP\P 

zuo 


9^ p.qr -4- n nnA 
zo.oyo in u.uu^ 


n 7kq _i_ n nno 

U. ( OO It U.UUO 


997 
ZZ / 


91 9 
Z1Z 


ZUl 


9^ ona -i- n nnA 
zo.yuo in u.uu^ 


n 787 4- n nnzt 

U.IOI It U.UU^t 


9n^ 

ZUO 


1 on 
iyu 


1 89 
loz 


op noo -i- n nri/i 
zo.uyy it u.uu4 


n 89fi 4 n nn^ 

U.oZO it U.UUO 


lyo 


1 8/i 

104 


1 79 
1 1 Z 


op 9ss -i- n nnA 

ZD. Zoo It U.UU^ 


n 84 <^ 4- n nnfi 

U.O4;0 It U.UUO 


1 7fi 
I/O 


100 


ioy 


9fi ^n9 -i- n nnzi 

ZO.OUZ It U.UU4: 


n so"? 4- n nna 
u.oyo it u.uuo 


ioy 


101 


1 ztzt 


op. pqr 4- n nn^ 
zo.uyo it u.uuo 


n ocifi 4 n nno, 
u.yoo it u.uuy 


1 37 
10/ 


loz 


1 OP 
1ZO 




c. HjI lQd-iiub 








9n ss7 -i- n n9K 

ZU.OO/ It U.UZO 


1 nn^ n 09^ 

l.UUO zn u.uzo 


/I 

4: 


A 
*1 


/I 


91 oon 4. n n99 
zi.ooy it u.uzz 


u.yoo zc u.uio 


■J 
O 


Q 
O 


Q 
O 


Zl.^O/ It U.UZ1 


u.yzo it u.ui4 











91 7m -i- n n9A 

Zl. / Ul It u.uz^ 


u.yio in: u.uuu 


S 



7 


7 


91 CQ9 _i_ n n^Q 
ZjL.ooZi it u.uoy 


U.o < O ZC U.Ulu 


O 


U 


O 


99 1 Qn _i_ n n9"? 

ZZ.lOU It u.uzo 


u.oyy in: u.uiz 


4: 


A 


/I 

4: 


99 9Szi -i- n n^A 

ZZ.Zo4 It U.UO^ 


u.o * o m u.uiu 


r. 
O 


r. 






99 ^9A -4- n m 8 

ZZ.OZ4 It U.Ulo 


U.oO ( HZ U.UUO 


Q 

y 








99 71 1 _|_ n m 8 
ZZ. / 11 It U.Ulo 


U.O^Z HZ U.UUO 


7 


(•. 



(•. 




99 RRP 4- n n99 
ZZ.OOV It u.uzz 


U.olo zt U.UUo 


1 9 
1Z 


1 1 
1 1 


1 1 
1 1 


9Q 1 n9 -i- n m n 

ZO.IUZ It U.U1U 


u.uyo zn u.uu * 


OO 


0^1 




23.303 ±0.009 


0.638 ±0.004 


40 


40 


40 


23.485 ±0.007 


0.624 ±0.004 


54 


53 


52 


23.710 ±0.007 


0.615 ±0.004 


59 


57 


54 


23.919 ±0.007 


0.621 ±0.004 


66 


64 


60 


24.094 ±0.007 


0.636 ±0.004 


82 


79 


76 


24.299 ±0.006 


0.646 ±0.004 


83 


80 


77 


24.496 ±0.006 


0.664 ±0.003 


125 


119 


111 


24.690 ±0.006 


0.682 ±0.004 


104 


97 


93 



-49- 



Table 2 — Continued 



V 


V-I 


ni 


^2 


n 3 


oa aaa -i- n nnfi 

Z4.ooo in U.UUO 


n 71 1 -i- n nnQ 

U. / 1 1 It U.UUO 


1 1 

110 


1 07 
1U 1 


1 m 

1U1 


zo.uyz it u.uuo 


n 701 _i_n 00^ 

U. / ol it U.UUO 


119 
1 1Z 


10^ 

1UO 


08 

yo 


ZO.OU/ in U.UUO 


n 7A7 4- nnfi 

U. ( 4 ( It u.uuo 


1 1 3 
110 


110 
1 1U 


1 07 
1U 1 


zo.ouz m u.uuo 


n 70,0 4- oo^. 
u. * yu it u.uuo 


1 

IOU 


1 O'i 
izo 


117 
1 1 1 


9c; 7m -i- n nn^ 
zo. / ui m u.uuo 


n si a -4- nnfi 

U.014 It U.UUO 


1 1 ^ 
110 


1 no 
iuy 


1 04 
1U4 


zo.oyy m u.uuo 


n szli 4- n nn7 

U.041 It U.UU ( 


1 41 
141 


1 LA 

104 


1 9Q 
izy 


op 1 na 4- o nn^ 
zo.iuo m u.uuo 


n orifi 4- n nn7 
u.yuo it u.uu * 


1 1 ^ 
110 


1 no 
iuy 


1 04 

1U4 


op. ^tnn -i- n nnfi 
zo.ouu m u.uuo 


n qao 4- n 000, 
u.y4z it u.uuy 


1 97 
1Z / 


1 91 
1Z1 


118 
llo 


op ^no 4- o nciP 
zo.ouz m u.uuo 


n cos 4- n nno 
u.yyo it u.uuy 


1 1 ^ 
110 


111 
111 


1 07 
1U 1 


op 7V\P 4- o nn7 
zo. / uo m u.uu / 


1 n*?^ 4- m 1 

1 .UoO It U.Ul 1 


Q1 

y 1 


8fi 
oO 


84 


zo.oyi m u.uuo 


1 ozls 4- ni n 

l.U4o It U.U1U 


8^ 

oO 


81 
ol 


an 

oU 


97 no9 -i- n nn7 
z * .uyz it u.uu / 


1 n/i 4- n m a 

l.U4o It U.U14 


OO 


PO 

oz 


oy 


97 97/L 4-d I1U 


1 no^ 4- n (19^. 
i.uyo it u.uzo 


9fl 

zu 


1 Q 

iy 


1 8 
lo 


97 zL79 -1- 097 
Z(.4/ZItU.UZ/ 


1 1 no 4- n ri9i 
i.iuy it u.uzi 


7 


7 


7 




H M 'i 
u. 1V1 








10.UO1 It U.UlO 


n QQ"? 4- O 004 

u.yyo it u.uu4 


Q 
O 


Q 
O 


Q 
O 


10.O4O It U.U14 


n ss9 4- oo^. 

U.OOZ It U.UUO 





O 


O 


1U.400 It U.Ul / 


n sal 4- oo^. 

U.ool It U.UUO 


Q 
O 


■3 



■3 



1 p pp(\ -4- n nn/i 

10.00U it U.UU4 


n 87/i 4- n ni 7 

U.o/4± U.Ul 1 


/I 


4 


4 


io.yoi it u.uzu 


n qap 4- n nna 

U.040 It U.UUO 


A 

A: 


A 

A: 


/I 


17.736 ±0.014 


0.810 ±0.007 


3 


3 


3 


17.931 ±0.010 


0.812 ±0.007 


5 


5 


5 


18.034 ±0.009 


0.799 ±0.007 


8 


8 


8 


18.235 ±0.007 


0.780 ±0.006 


8 


7 


7 


18.363 ±0.007 


0.743 ±0.008 


19 


18 


18 


18.455 ±0.006 


0.680 ±0.009 


20 


20 


20 


18.552 ±0.007 


0.647 ±0.005 


23 


22 


21 


18.655 ±0.009 


0.608 ±0.011 


14 


14 


14 
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V 


V-I 


m 


ri2 


n 3 


lo. / 44 in U.UUO 


n fiin + n nn"? 
u.oiu m u.uuo 


3A 
04 


39 
oz 


3n 
ou 


1 s sr^ -i- n nn^ 

lo.oOO in U.UUO 


n KQ'i 4- n nn"? 
u.ooo m u.uuo 


^n 
ou 


9Q 
zy 


9Q 
zy 


i s q^"? -i- n nnR 
lo.yoo in u.uuo 


n krk 4- n nn9 
u.ooo m u.uuz 


A\ 


3Q 

oy 


37 
t 


1 q nA8 -i- n nnA 
iy.u4o m u.uu4 


n -i- n nn9 
u.ooi m u.uuz 


AA 
44 


A1 
41 


An 

4U 


iy . 101 m u.uuo 


n kra -4- n nn9 
u.oo4 m u.uuz 


GO 


ou 


AQ 
4y 


1 q 9^n -i- n nn"? 
ly.zou m u.uuo 


n -i- n nn9 
u.ooi m u.uuz 


^.Q 

oy 


^A 
04 


01 


1 q "?ar -+- n nnA 
iy.040 m u.uu4 


n ^si -i- n nn9 
u.ooi m u.uuz 


00 


ou 


A^ 
40 


1 q AA7 4- n nnA 


n kqq 4- n nn9 

u.ooy m u.uuz 


oo 


AQ 
4y 


AQ 

4y 


1 q ^a^ -i- n nnA 
iy.040 m u.uu4 


n c;on -i- n nn9 
u.oyu m u.uuz 


Rn 
ou 


^7 
O ( 


^A 
04 


1 q r^r -i- n nnA 
ly.uoo m u.uu4 


n gar -i- n nn^t 
u.uuu m u.uuo 


oo 


^9 
OZ 


^n 
ou 


1 q 7^9 -i- n nnA 
iy. / oz m u.uu4 


A GAQ _l_ A AAQ 

u.uuy m u.uuo 


00 


oo 


AQ 
4y 


1 q sas -i- n nn^t 
iy.o4o m u.uuo 


A fil A _l_ A AAO 

u.014 m u.uuz 


RR 
OO 


R 1 
1 


^R 
oo 


1 q q^a -i- n nnA 
iy.yo4 m u.uu4 


A «Ofi _|_ A AAQ 

u.uzo m u.uuo 


R7 
O 1 


R'i 
OO 


Rl 
01 


9n n^7 -i- n nnA 
zu.uo/ it u.uu4 


A «QQ _|_ A AAO 

u.ooo m u.uuz 


73 
i 


RQ 

oy 


R^ 
OO 


on i aq -i- n nnA 
zu.i4y m u.uu4 


n p,Ap, 4- n nn"? 
u.040 m u.uuo 


7A 
( 4 


RQ 

oy 


R3 
oo 


9n 9A7 -i- n nni 
zu.z4/ m u.uuo 


A GEQ _l_ A AAO 

u.ooy m u.uuz 


sn 

OU 


7 s . 
( 


79 
/ z 


9n iaq -i- n nni 
zu.04y m u.uuo 


A G70 _l_ A AAO 

u.o / z m u.uuo 


78 


77 
i i 


7A 
/ 4 


9n a^a -+- n nni 
zu.404 m u.uuo 


n «Q7 _i_ a nn9 
u.oo / m u.uuz 


89 
oZ 


7^ 
1 


7n 


9n ka r -i- n nn/i 

ZU.040 it U.UU4 


n 7ao _i_ a nn°. 

u. / UZ it U.UUo 


Q7 

y t 


an 
yu 


8R 
oO 


9n rar -i- n nni 
zu.040 m u.uuo 


n 7it; + n nni 

U. / 10 It U.UUO 


Q3 

yo 


Q1 

y i 


8A 
o4 


20.754 ±0.004 


0.734 ±0.004 


89 


88 


82 


20.850 ±0.004 


0.749 ±0.003 


64 


59 


55 


20.950 ±0.003 


0.760 ±0.003 


95 


90 


84 


21.053 ±0.004 


0.777 ±0.003 


85 


81 


77 


21.149 ±0.003 


0.796 ±0.003 


75 


72 


67 


21.247 ±0.003 


0.813 ±0.005 


67 


65 


61 


21.352 ±0.004 


0.829 ±0.003 


58 


55 


53 


21.444 ±0.004 


0.847 ±0.004 


61 


56 


52 
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V 


V-I 


m 


^2 


«3 


Zl.OOZ In u.uuo 


n R7K _i_ n nn"? 
u.o / m u.uuo 


ZLQ 
4y 


zlfi 

40 


zLzL 
44 


zi.04y it u.uu4 
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u.oyo it u.uu4 


ZlQ 

4y 


A K 
40 
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4Z 


oi _i_ n nn^ 

Zl. / 00 in U.UUO 


n 01 a 4- n nnzL 
u.y 14 m u.uu4 


OO 


ou 


ZL7 
4 / 


oi pin _i_ n nnzL 
zi.o4y in u.uu4 


n cm 1 -i- n nn^ 
u.y4i m u.uuo 


ZLQ 
4y 


ZL7 
4 / 


40 


91 007 _i_ n nn^ 
zi.yo/ m u.uuo 


u.yoo m u.uuo 


44 


ZL1 
41 


oy 


zz.uoo in U.UU4 


n 0Q7 _i_ a nnK 
u.yo / m u.uuo 


ZLQ 
4y 


ZL7 
4 1 


zt^ 

40 


99 uoxn nn^ 
zz.i^z in u.uuo 


1 nno -i- n nn^ 
i.uuy m u.uuo 


ZL1 
41 


00 


OO 


99 oka _i_ n aac 

zz.zj^i in u.uuo 


1 aqq _i_ a nnfi 

i.uoo in u.uuo 


OO 


"?9 
oz 


01 


99 'ZA'z -4- n nn^ 
zz.o^o m u.uuo 


1 n^ -i- n nn^ 
i.uoo m u.uuo 


oy 


OO 


^ZL 
04 


99 zL^n -i- n nnzL 
zz. 40U m U.UU4 


1 n7s _i_ a nnzL 
i.u / m u.uu4 


ztfi 

40 


49 
4Z 


ZL1 
41 


99 kkci 4- n nn^ 
zz. oou m u.uuo 


1 1 1 Q _l_ A A AO 

1.110 in u.uuo 


01 


90 
zy 


97 
Z / 


99 rv7 -4- n nns 
zz.oo/ in u.uuo 


1 1 8Q _|_ A A1 K. 

1 . loo in u.uio 


1 7 
1 1 


10 


1 f; 
10 


99 -i- n nno 
zz. / 00 m u.uuy 


1 1 ftft _l_ A A9I 

1.100 in u.uzi 


1 n 

1U 


1 n 

1U 


i n 

1U 


99 s^n -i- n nno 
zz.oou m u.uuy 


1 1 7s 4- n m q 
1.1 / in u.uiy 


a 











C. 1V1 








1 /i qka 4- n n9i 
i4.yo4 m u.uzi 


1 nni -4- n m 9 
1.UU4 m u.uiz 


Q 
O 


q 



Q 
O 


lo.uoy m u.uoo 


A QQG _l_ A AAO 

u.yyo m u.uuo 


q 
O 


O 


Q 
O 


10.0 / u m u.uio 


A CM fi + n 111"! 

u.y4o m u.uio 


/I 
4 


A 
4 


A 
4 


10.040 it U.UoZ 


A QA C _l_ A AAA 

u.y4o in u.uuy 



O 


q 



q 



1 c; qka _i_ n n9i 
10.004 in u.uzi 


A QO A _l_ A A1 1 

u.you m u.ui 1 


/I 
4 


4 


A 
4 


15.965 ±0.016 


0.928 ±0.011 


4 


4 


4 


16.432 ±0.013 


0.894±0.012 


5 


5 


5 


16.557±0.017 


0.890 ±0.006 


5 


5 


5 


16.678 ±0.007 


0.888 ±0.005 


4 


4 


4 


16.758 ±0.005 


0.873 ±0.005 


5 


5 


5 


16.948 ±0.009 


0.877 ±0.008 


8 


8 


8 


17.042 ±0.013 


0.869 ±0.010 


7 


7 


7 


17.156 ±0.005 


0.865 ±0.004 


12 


12 


12 
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V 


V-I 


m 


^2 


n 3 


1 7 98"? -i- o oo^ 
i / .zoo m u.uuo 


n 814 -i- ooq 
u.ooi m u.uuy 


Q 
O 


Q 
O 


Q 
O 


1 / .001 in U.U10 


n 848 -1-0 011 
u.o4o in u.uu 


fi 



fi 
O 


fi 



1 7 419 -i- n m i 

± / .40Z in U.U11 


A 8^8 -1- 008 

u.ooo in u.uuo 


K 
O 


K 
O 


r. 
O 


i 7 rcrco -1-0 010 
1 / .000 in U.U1U 


n 898 -1- 007 
u.ozo in u.uu/ 


1 
1U 


1 n 

1U 


1 
1U 


1 7 K^l -4- 000 

± / .000 m u.uuy 


n 830 -1- 007 
u.oou in u.uu / 


1 1 

1 1 


1 1 

1 1 


1 1 

1 1 


1 7 7CI -1-0 011 


n soQ -4-0 01*? 
u.ouy m u.uio 


fi 



7 


7 


1 7 S^fi -1- 007 

± / .000 m u.uu / 


n 7R7 -I- OOfi 

u. / / in u.uuo 


9fi 
ZO 


9fi 
ZO 


9fi 
ZO 


1 7 qko _i_ n oofi 
1 f .yoz m u.uuo 


n 7i74-n 00^ 

u. / 1 / m u.uuo 


14 
04 


11 
Ol 


10 
OU 


1 8 0^0 -I- OOfi 
lo.uou in u.uuo 


n R74 4- 00"? 
u.o / 4 m u.uuo 


97 
Z / 


9 s . 

zo 


91 

zo 


18 1 48 -1- 007 
±0. 140 m u.uu / 


A f{A 1 _i_ 009 
u.041 m u.uuz 


in 
ou 


97 
z / 


9fi 

zo 


1 8 948 -1- OOR 
io.z4o m u.uuo 


a coo -1- 001 
u.oou m u.uuo 


ifi 
00 


19 

oz 


11 

Ol 


1 e q«i _i_ OOzl 
±0.001 m U.UU4 


A f?94 -1- 00"? 

u.oz4 m u.uuo 


18 
00 


17 
/ 


IK 
OO 


1 8 4^1 -I- 00^ 
10.400 m u.uuo 


A GOI _i_ 003 

u.001 m u.uuo 


4fi 
40 


4 s . 

4:0 


49 
4Z 


1 8 ^17 -I- OOfi 
10.00/ m u.uuo 


a «o/i _i_ 00 s ; 
U.OZ4 m u.uuo 


9fi 

zo 


9 s . 

zo 


94 
Z4 


1 8 RA<^ -4- OOzL 
10.040 m U.UU4 


A G97 _i_ 009 
u.oz / m u.uuz 


^.8 

DO 


OO 


^1 
OO 


1 8 7/1Q -1- 004 

10. / 4y m u.uu4 


fill -1- 001 

u.001 m u.uuo 


4Q 
4y 


4fi 


4^ 
40 


1 8 847 -1- 004 


fill _l_ A AAO 

u.ooo m u.uuo 


^0 

ou 


40 
4y 


4fi 
40 


1 8 0^8 -1- 00^ 

lo.yoo m u.uuo 


fi41 -I- 001 
u.041 m u.uuo 


48 
4o 


4fi 
40 


44 
44 


1 q -i- nm 
ly.UOO it u.uuo 


R48 -1- 001 

u.04o in u.uuo 


^O 
OU 


/1Q 
4y 


1 8 

4o 


in uoj_n 00^ 
iy.100 m u.uuo 


p,KC\ -4- 009 

u.oou m u.uuz 


40. 

4y 


48 
lo 


4^ 
40 


19.250 ±0.004 


0.661 ±0.003 


56 


55 


53 


19.352 ±0.003 


0.668 ±0.003 


74 


70 


69 


19.450 ±0.004 


0.675 ±0.002 


60 


58 


54 


19.551 ±0.004 


0.688 ±0.002 


89 


86 


82 


19.658 ±0.004 


0.698 ±0.003 


50 


48 


47 


19.751 ±0.004 


0.709 ±0.003 


69 


67 


64 


19.851 ±0.004 


0.719 ±0.003 


65 


62 


61 


19.955 ±0.004 


0.742 ±0.003 


60 


56 


54 
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V V-I n\ ri2 nz 



20.049 ±0.004 
20.154±0.004 
20.255 ±0.004 

20.345 ±0.004 
20.448 ±0.003 
20.545 ±0.004 
20.649 ±0.004 
20.744 ±0.004 
20.852 ±0.002 
20.946 ±0.004 
21.044 ±0.005 
21.140 ±0.004 
21.258 ±0.005 
21.356 ±0.005 
21.441 ±0.004 
21.552 ±0.005 
21.660 ±0.005 
21.746 ±0.006 
21.858 ±0.007 
21.958 ±0.006 
22.058 ±0.007 
22.134 ±0.007 
22.251 ±0.011 

22.346 ±0.015 
22.537 ±0.023 



0.746 ±0.003 
0.765 ±0.003 
0.780 ±0.003 
0.783 ±0.004 
0.809 ±0.003 
0.829 ±0.004 
0.840 ±0.003 
0.857 ±0.005 
0.875 ±0.003 
0.898 ±0.005 
0.916 ±0.004 
0.930 ±0.004 
0.969 ±0.005 
0.985 ±0.005 
1.002 ±0.004 
1.036 ±0.006 
1.058 ±0.006 
1.073 ±0.010 
1.095 ±0.009 
1.162 ±0.007 
1.168 ±0.010 
1.186 ±0.013 
1.275 ±0.017 
1.244 ±0.043 
1.350 ±0.024 



85 


81 


77 


65 


63 


60 


63 


62 


60 


58 


56 


52 


66 


63 


62 


62 


59 


57 


64 


60 


58 


60 


57 
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Table 3. Turnoff Locations for Program Clusters 



Cluster Vmsto (V-I)msto V+o.io V+om ^+0.05 V+0.10 



Palomar 3 


23.83 


0.614 


23.22 


23.35 


24.63 


25.09 


Palomar 4 


24.12 


0.614 


23.52 


23.64 


24.93 


25.38 


Eridanus 


23.70 


0.610 


23.09 


23.21 


24.50 


24.95 


M3 


19.16 


0.578 


18.53 


18.67 


19.97 


20.40 


M5 


18.56 


0.614 


17.95 


18.09 


19.37 


19.84 
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Table 4. Giant-Branch Offsets for Program Clusters 



mean median standard median 



Cluster 


offset 


offset 


deviation 


N 


a( V-I) 


Palomar 3 


-0.0035 


-0.001 


0.021 


45 


0.013 


Palomar 4 


-0.0037 


-0.003 


0.012 


51 


0.011 


Eridanus 


+0.0078 


+0.009 


0.020 


46 


0.014 
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Table 5. Giant-Branch Offsets for Program Clusters 



mean median standard median 



Cluster 


offset 


offset 


deviation 


N 


a( V-I) 


Palomar 3 


+0.0096 


+0.015 


0.020 


44 


0.013 


Palomar 4 


+0.0123 


+0.012 


0.010 


50 


0.011 


Eridanus 


+0.0205 


+0.026 


0.022 


45 


0.014 


M3 


-0.0139 


-0.015 


0.014 


46 


0.005 


M5 


-0.0073 


-0.006 


0.019 


85 


0.010 
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Table 6. Fiducial Parameters for Program Clusters 



Cluster 


[Fe/H] 


E(B-V) 


(V-/) ,to 


V HB 


Vmsto 


AV 


Palomar 3 


-1.57 


0.03 


0.579 


20.51 


23.83 


3.32 


Palomar 4 


-1.28 


0.01 


0.595 


20.80 


24.12 


3.32 


Eridanus 


-1.42 


0.02 


0.590 


20.42 


23.70 


3.28 


M3 


-1.57 


0.01 


0.578 


15.68 


19.14 


3.46 


M5 


-1.29 


0.03 


0.579 


15.15 


18.59 


3.44 




Stetson et al. Fig. 2 




Stetson et al. Fig. 5 
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Stetson el al. Fig. 13 
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Stetson el al. Fig. 15 
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Fig. 19 
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